Exploring initial correlations in a Gibbs state by local application of external field 
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We demonstrate that local application of a weak external field increases distinguishability between 
states with and without initial correlations. We consider the case where a two-level system linearly 
and adiabatically interacts with an infinite number of bosons. We evaluate the trace distance (or, 
equivalently, the Hilbert-Schmidt distance) between the quantum states which evolve from two kinds 
of initial states; the correlated Gibbs state and its uncorrelated marginal state. We find that the 
trace distance increases above its initial value for any and all parameter settings. This indicates 
that we can explore the existence of initial system-environment correlations with applying a local 
external field, which causes the breakdown of the contractivity. 

PACS numbers: 03.65.Yz, 03.65.Ta, 03.67.-a 



I. INTRODUCTION 

Reduced dynamics of open systems plays an important 
role in various fields such as condensed matter, quantum 
optics, and quantum electrodynamics In order to de- 
scribe the reduced dynamics, an initial condition is often 
used where the relevant system is statistically indepen- 
dent from its environment. The uncorrelated condition 
was originally used to describe NMR (nuclear magnetic 
resonance) phenomena for systems with weak system- 
environment interactions 0, Besides the initial con- 
dition, the Markovian approximation has been used to 
describe the relaxation phenomena, which gives the so- 
called Bloch equation^]. Completely positive dynamical 
semigroups have been devised to ascertain decay proper- 
ties in the Markovian process [5], imparting contractivity 
to the open system dynamics. While dynamical semi- 
groups became the starting point for design of quantum 
information processing^], as experimental techniques de- 
veloped to observe non-Markovian dynamics [7[, theoret- 
ical treatments for such behavior [l| attracted much at- 
tention. The projection operator methods have provided 
formulations for initially correlated systems [a-tlQj]. but 
the initial correlations have been frequently ignored (Tl). 

The effects of initial system-environment cor- 
relations have been discussed regarding quantum 
measurements II 2h24| . complete positivity of dynam- 
ical maps [25Tt3l| . non-Markovian dynam ics [ 33 -l34| , 
linear-response absorption line shapes [351 - 1371 ] . and 
general theoretical treatments using the projection 
operator methods (38j. Especially relevant has been 
the finding that initial system-environment correlations 
require extension of the conceptual framework of open 
system reduced dynamics to include generalization of 
contractivity (39j . 

The breakdown of contractivity by including ini- 
tial correlations has been examined by evaluating the 
time evolution of the trace distance for finite size of 
environments [39l l4fj| . which shows an increase over the 
initial values. Dajka, Luczka, and Hanggi compared 
four different distance measures: the trace (which is 



equivalent to Hilbert-Schmidt), Bures, Hellinger, and 
quantum Jensen-Shannon distances, and found that only 
the trace distance reveals the distance increase for a 
system that initially correlates with an infinite size of 
environment [4l| . However, the correlated state in [4l[ 
requires elaborate preparation by quantum engineering. 
As a correlated state, we can consider a Gibbs state 
which includes infinite size of environment. When we 
analyze NMR and ESR (electron spin resonance) experi- 
ments, or design quantum information processing for con- 
densed matter, it might be necessary to find the system- 
environment correlations that are inherent to the mat- 
ter. A recent study on the transient linear response of a 
matter system to a suddenly- applied weak external field 
showed dependence on the initial conditions, correlated 
Gibbs state, and a conventional factorized state where 
the relevant system and the environment stay in their 
equilibrium states [42{. Conversely, one might say that 
the transient linear response of a matter system is use- 
ful to detect the initial correlations. However, in [42j . 
it was difficult to distinguish the effect of initial correla- 
tions except for strong system-environment interactions 
and for intermediate temperatures. In addition, within 
this approach, we cannot discuss the generalization of 
contractivity. 

In this paper, we evaluate the time evolution of trace 
distance for a matter system where a two-level system 
interacts with an infinite number of bosons as the en- 
vironmental system. We consider that the two-level sys- 
tem linearly and adiabatically interacts with the environ- 
ment, and that a weak external field is suddenly applied 
to the two-level system. We show that the trace dis- 
tance between the time evolution from the Gibbs state, 
Pse, and from an uncorrelated initial condition clearly 
exhibits effects of initial correlations as an increase in 
the short time region. As opposed to [42{ . as an uncor- 
related initial condition, we consider the marginal state 
of the Gibbs state, namely ps ® Pe with ps — T^ePse 
and pe — TtsPSE where Trs(TrE) is a partial trace over 
the relevant (environmental) system, respectively. By de- 
scribing an uncorrelated state as the marginal state, the 
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relevant system interacts with the same environmental 
features as an average. We find that the trace distance 
increases in the short time region even when the system- 
environment interaction is weak. This indicates that the 
local application of a weak external field can cause in- 
formation inaccessible at an initial time to flow into the 
relevant system from its infinitely-sized environment in- 
cluded in a Gibbs state, and we can monitor it with the 
trace distance effectively. 

The paper is organized as follows. In Sec. II, we pro- 
vide our formulation to obtain the time evolution of the 
two-level system which linearly and adiabatically inter- 
acts with an infinite number of bosons under a suddenly- 
applied weak external field. We evaluate the induced 
dipole moment of the two-level system in Sec. Ill, and 
the trace distance time evolution from two different ini- 
tial conditions in Sec. IV. We state conclusions in Sec. 
V. 



II. FORMULATION 

We consider a matter system which consists of a two- 
level system and its environment of an infinite number of 
bosons. We assume that the system-environment inter- 
action is linear and adiabatic, which causes the pure de- 
phasing phenomena to the two-level system. The Hamil- 
tonan of the matter system is written as 



the matter Hamiltonian in the form 

W = S^HS = H' S + H E , (4) 
H's = E[\1){1\+E \0)(0\, (5) 

with E[ = Ei — fij^kidk/^k)- The matter-field interac- 
tion is transformed as 

U' P {t) = S^H P {t)S 

= --p-H9(t)\l)(0\e- iu ''' t G{Q) + h.c., (6) 

with G(0) = e s . From Eqs. (j4|)- ((6]), the time evolution 
of the matter system takes the form 

P (t) = U{t) P {0)U\t) = S , e-* w 'V(t)^(0)Vt(t) e *w'*St j 

where we define the transformed initial condition as p'(0) 
and 

U{t) ee T + exp[-{^Ht + J dt' l -U P {t'))l (8) 

V(t) ee T + cxp[-^J*<H' P (t>)dt'} (9) 

with H' p {t) ee e*' H ' t n p (t)e--k' H ' t {4:3}. In the following, 
we consider the application of a weak external field, which 
enables us to approximate, up to first order, the matter- 
field Hamiltonian in V(t) as V(t) w 1- | J*H' P (t')dt' . 



with 



H = Hs + He + H S e, 



Us = £i|l)(l|+£o|0><0|, 

7-L E = fjLdkbjbk, 
k 

n SE = + 



(i) 



(2) 



where Eq and E\ are the energy of the lower and up- 
per states, respectively, of the relevant system, Uk is the 
frequency of the bosonic bath mode, bt and bk are its cre- 
ation and annihilation operators, and gk is the coupling 
strength. We study the transient behavior of a two-level 
system after a sudden application of an external field. 
The matter-field interaction Hamiltonian is given by 



1 



H P (t) = --p-H0(t)\l)(O\e 



1 



where p, is the transition dipole moment, El is the am- 
plitude of the external field, uj p is the frequency of the 
external field, and 9(t) is the step function. 

In order to evaluate the time evolution of the induced 
dipole moment, we use a canonical transformation in 

terms of 5 ee exp[B|l)(l|] with B = Efc(.9fcM)( & fc " b l), 
since we can eliminate the system-bath interaction from 



A. correlated initial condition - Gibbs state - 

As a correlated initial condition, we consider that 
the whole matter system is in a Gibbs state, p(0) = 
Pse = h exp[— f3TL\. Transforming the correlated ini- 
tial state, we find p'(0) = tf cxp[— PU' S \ exp[— (3He] with 
Z = Tr[exp[— j3"H]] where Tr is the trace operation for 
the total system, which gives 



p , (0)=S^(0)S = p' o \0){0\+p' 1 \l){l\ 



with 



Pi) 



Pi = 



e-^PE 



Z's 



(10) 



(11) 



In Eq. (HD, we define Z' s = Tr s [exp[-/37^]] and p E = 
exp[-(3n E }/Z E with Z E = Tr E exp[-pH E ]. Note that 
the transformation of the correlated state gives a fac- 
torized state of the Gibbs state and the bosonic environ- 
ment, and endows the two-level system with renormalized 
energies. 

The transformed initial condition gives the elements of 
the reduced statistical operator as 



Tr E [U(t)p SE U(t)} 
TiE[U(t)p SE U(t)} 



-f3E[ 



7' ' 

ip, ■ H 



2h 



C 4(corr)(t),(12) 
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where we define 



A (corr) (i) = / df e -"M*-0 



X {-=r-*i(t t') -=-**(* t% (13) 

with Aw = (E[ - Eq)/% - U) p . In Eq. (fl3]). we define 

*i(t - = (Gt(t)G(i')) = (G(t')Gt (<))*, (14) 

with G(t) = e^e^V^'*, and (X) = Tr E/ 5 B X for an 
arbitrary operator X. 

B. uncorrelated initial condition - marginal state - 

As an uncorrelated initial state, we take a factorized 
state, p(0) — ps ® Re which is the marginal of pse with 
ps = T?ePse and pE — Ti'sPSE- (Note the difference 
between pe, and p~e)- As shown in the Appendix, the 
transformation of the marginal state gives 

Po = ^re-P E °p E , p[ = ie-^G(0)p B Gt(0), (15) 

with p E = ^r(e-P Eo p E + e-^Gt(0)/5 E G(0)). The el- 
ements of the reduced statistical operator are obtained 



Tr E [U(t) Ps ®p E U(t)} 
TMU(t)ps®p E U(t)} 



-f3E[ 
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'j[_lL, .... 

2R 



^(marg) (*)) 

(16) 



where we define 



A 



(marg) 



(*) 



d p e -iAQ(t-t') 



p -0E o p -j3E _ p-PK - 



using the formulation in the previous section. Apart from 
a factor of 2p(jj,* ■ H)/H, we define the induced dipole 
moment for the correlated initial condition as /X( corr ) (t) 
and for the uncorrelated marginal initial condition as 
M(mar g )(*)- We obtain 

= \ A (m)(t)\ COs(uj p t - </>( m )(f)), (20) 



where A( corr )(t) and Ar maTS )(t) are defined in Eqs. (fT3|) 
and (|17[) . respectively. (f)/ m \(t) is the argument of 
A( m )(t). Defining the coupling spectral density as h(uj) = 
J2k 5^0^ _ ^k), we obtain 



and 



exp[ 



da; 



x{(l + 2n(w))(l - cos(urf)) +isin(wt)}], 

(21) 



* 2 (t,t / ) = 

x exp[— / da; — — 2z(sin(u;t') - sin(a;t))l. 
Jo w 

(22) 

Let us now set the spectral density to be Ohmic, i.e., 



h(u) 



j/u 



: , with the coupling strength, 5, and 



(17) 



with 



¥ 3 (t,f) ee (G(0)Gt(t)G(t')Gt(0)) 
= (G t (0)G(t')G t (i)G(0))* 



(18) 



III. TIME EVOLUTION OF INDUCED DIPOLE 
MOMENT 

We evaluate the induced dipole moment under the ap- 
plication of an external field given by 



/2(t) = Tr[(/2|l)(0| + /r|0)(l|)p(t)] 



(19) 



the cut-off frequency, u) c . In this case, the renormalized 
energy is given by E[ = E\ — Hsuj c , and we find 

\t 2 (t,i') = &i(t—t') exp[-2is(ajctan(w c f')-arctan(a; c t))]. 

(23) 

In Fig.[T] we show the time evolution of the intensity of 
the dipole moment under the application of an external 
field for kgT = IOTiloq, s = 1, u) c = 0Jq/5, and lo p = luq. 
Time is scaled as t = uoqI. Dashed (blue), and light 
gray solid (orange) lines represent the time evolutions of 
induced dipole moment for the correlated and marginal 
initial conditions, respectively. 

We find that the induced dipole moment for the corre- 
lated initial condition approaches a stationary oscillation 
faster than for the marginal initial condition. As shown 
in Fig. QTb) and (c), the phase and amplitude of the in- 
duced dipole moment for the marginal initial condition 
approaches the same phase and amplitude of the cor- 
related initial condition. We can explain the long-time 
behavior using the monotonic time dependence of the 
arctan (u> c t) function. The function approaches 7r/2 with 
increasing time, which means that "^(M') approaches 
&i(t — t') for large t. Comparison between Eqs. (fT3|) 
and (|17[) reveals that the induced dipole moment for the 
correlated initial condition agrees with that for marginal 
initial condition at long times. 

The time evolutions of the induced dipole moment fol- 
lower temperatures of ksT = Hloq and fc^T ' = huio/5 are 
shown in Figs. [2] and [3l respectively. As the tempera- 
ture decreases, we find that the induced dipole moment 
for the correlated initial condition approaches a station- 
ary oscillation slower, and find a decreased dependence of 
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FIG. 1: (color online) Time evolution of the induced dipole 
moment for ksT — 10/iajo, s = 1, u} c = (jJo/5, and uj p = ujo- 
Dashed (blue), and light gray solid (orange) lines repre- 
sent the time evolutions of induced dipole moment for cor- 
related and marginal initial condition, respectively. Time 
is scaled as t — uiot. (a) time evolution of intensity of 
dipole moment |/i( m )(t)| 2 , (b) the time evolution of ampli- 
tude |A( m )(i)| and (c) the time evolution of phase, cf>r m \{t) 
with m = corr or marg. 

the time evolution on the initial condition. In the above 
evaluation, the stationary oscillation for the marginal ini- 
tial condition coincides with that for the correlated initial 
condition. This coincidence is in contrast to the behav- 
ior seen with the factorized initial condition, where the 
system and environment each stay in equilibrium, as dis- 
cussed in |42f . 

IV. TRACE DISTANCE 

We found the difference between the time evolution of 
the induced dipole moment from the Gibbs state of the 
whole matter system and its marginal state in the previ- 
ous section. However, it might be difficult to observe the 
difference between these time evolutions clearly. More- 
over, we cannot discuss contractivity from the evalua- 
tions in the previous section. 

In order to clarify the effects of initial system- 
environment correlations on open system dynamics, we 
need a sensitive and tractable measure. Since the ef- 
fects of initial correlations appear in time evolutions 
where non-Markovian features dominate, we can use 
the trace distance as a measure of non-Markovianity 
0, HH . The trace distance for two quantum states ex- 
pressed with trace class operators p 1 and p 2 is defined as 
D{p 1 1 p 2 ) = ^Trlp 1 — p 2 |[6j. Here we want to obtain the 
distance between the reduced dynamics of the two-level 
system for two kinds of initial conditions: the Gibbs state 



FIG. 2: (color online) Time evolution of the induced dipole 
moment for fcsT = Tiojq. Other parameters and evaluated 
quantities are the same as in Fig. [T] 
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FIG. 3: (color online) Time evolution of the induced dipole 
moment for fcsT = huuo/5. Other parameters and evaluated 
quantities are the same as in Fig. [T] 



p 1 (0) = pse, as a correlated initial condition, and its 
marginal state, p 2 (0) = ps®PE, as an uncorrelated initial 
condition. When the trace distance increases above an 
initial value, namely, D{p 1 {t), p 2 (tj) > D{p 1 {Q), p 2 (0)), 
we deduce that information, which is inaccessible at 
the initial time, flows into the relevant system through 
system-environment correlations. Moreover, this increase 
indicates a breakdown of the contractivity [39j. 
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For our model, we obtain the trace distance in the form 

D{Tr E [U(t) PSE U f (t)] , Tr E [U(t)p s ® (t)] ) 

= v/IXnP + l^iol 2 , (24) 

with X = Tifi^CtJpsB^W] ^ Tr E [i7(t)p5 ® pst^(t)]. 

Using Eqs. (|12p and Eq. (|16l) . we obtain the trace dis- 
tance as 

DCTrsIl/tpsBl/j^TrK^tPflriapB^]) - |Mio|- (25) 

We show the time evolution of the trace distance, apart 
from the dimensionless quantity J^—, for Ohmic spectral 
density in Fig.|4]by setting uj c = wo/5, and co p = loq. Here 
we scaled time variable as t (= uiot). 

Figures HJa)-(c) correspond to the trace distance time 
evolution at various temperatures: (a)fcsT = IOHujo, 
(b)k E T — Hujq, and (c)k E T = Tuo^/h. In each figure, 
we vary s as 1, 0.1, and 0.05. We find in Fig. QJa) 
that the trace distance increases in the short time re- 
gion, which signifies the breakdown of contractivity for 
each value of s. In Fig. Ufa), we also find the trace dis- 
tances approach zero at long times, which arises from the 
fact that the arctan (w c t) function in Eq. (|23p monoton- 
ically approaches 7r/2 with increasing time and ^(i, i') 
approaches ^i(t — t') for large t. For lower temperature 
cases as shown in Figs.[4|b) and (c), we find qualitatively 
similar behavior as in Fig.|3Ja). Comparing these figures, 
we find that the peak value of the trace distance becomes 
larger for increasing values of s and at intermediate tem- 
peratures. 

In the above evaluations, we find that the reduced dy- 
namics of the matter system under a weak external field 
shows a breakdown of the contractivity for any param- 
eter setting, which allows us a method to monitor the 
initial system-environment correlations in a Gibbs state. 
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FIG. 4: (color online) Time evolution of the trace distance 
for cj c = Wo/5, and cu p — luq with changing temperature as 
(a)fcsT = lOftcjo, (b)fcsT = huio, and (c)fcsT = ?iajo/5. In 
each figure, black, light gray (red), and dark gray (blue) lines 
correspond to the case of s = 1, 0.1, and 0.05, respectively. 
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V. CONCLUSIONS 

We have studied the time evolution of a matter system 
which is composed of a two- level system and a bosonic en- 
vironment of infinite size, under a suddenly-applied weak 
external field. Assuming the system-environment inter- 
action to be linear and adiabatic, we evaluated the in- 
duced dipole moment for two kinds of initial conditions: 
a correlated Gibbs state and an uncorrelated marginal 
state. We found that the time evolution of the induced 
dipole moment depends on the initial condition. How- 
ever, the difference depends on the parameter settings 
and, under certain settings, may be difficult to observe. 
By evaluating the trace distance in the reduced dynam- 
ics for these two initial conditions, we found that we 
can overcome the parameter dependence. The trace dis- 
tance clearly increases in the short time region for cases 
of arbitrary system-environment interactions and tem- 
peratures. We can effectively monitor the initial system- 
environment correlations in a Gibbs state with the break- 
down of the contractivity. 



Appendix: Derivation of Eq. (|15[) 

In this appendix, we explain how to obtain the 
marginal initial state 

p{0)=p s ®PE (AT) 

with p s = TtePSE and p E = Trsp.SE for pse{= 
A exp[— &%}). The correlated state pse is rewritten as 

Pse 

1 r p 

= - exp[-/3^ ]T+ exp[- / d\e XHo Hs E e- XH "] 
^ Jo 

= I(e-^o exp[ _ W | 0)(0 | 
+e- 0El exp[-(3H E ] 

xT + cxp[- V% fe [ d\(ble hXu * + b k e- hXu *)]\l)(l\), 
k Jo 

(A.2) 
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where we define Hq = Tis + T~Le and use the T + operator From these, we obtain the marginal state as 
to order A from right to left. Using the relation as 



exp[-/3H B ] T+ exp[- V fig k f d\{b\e hx " k + b k e 
k Jo 

= exp[-/3(H B + n 9k{b\ + b k ))} 



PS ® PE = ^-(e-"*>p B |0)<0| + e-^p B |l)(l|), (A.7) 



G*(0)exp[-p(H E -hY^)]G(0), 



(A.3) 



we obtain 



Pse = -^(e^ E «p E \0)(0\ + e-^Gt(0)p E G(0)|l)(l|), 

(A.4) 

which gives the partial trace operation on p$E over the 
system and the environment as 

PS = i(e-^|0)(0| + e-^|l)(l|), (A.5) 



which is transformed into 



(ps®Pe)' = ^|0)(0|+pi|l)(l| 

= W !iEa PE 
z s 



+e-^G(0)p E Gt(0)|l)<l|). (A.8) 



p E = 4r(e~ PEo PE + e-^ B; G t (0) / 5 B G(0)). (A.6) Equation (|X8]) gives Eq. (JT5 
Z s 
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